REGULARIZATION OF THE KEPLER PROBLEM ON THE SPHERE 



SHENGDA HU AND MANUELE SANTOPRETE 



Abstract. In this paper we regularize the Kepler problem on S 3 in several different ways. 
First, we perform a Moser-type regularization. Then, we adapt the Ligon-Schaaf regularization 
to our problem. Finally, we show that the Moscr regularization and the Ligon-Schaaf map 
we obtained can be understood as the composition of the corresponding maps for the Kepler 
. problem in Euclidean space and the gnomonic transformation. 
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1. Introduction 



It often happens that the flow associated to a vector field is incomplete. Famous examples 
are the Newtonian n-body problem, where singularities arise because of the existence of colli- 
. sion orbits, and the Kepler problem. Regularization of vector fields is a common procedure in 
the study of differential equations. There are two main approaches. In the first approach the 
incompleteness of the flow is removed by embedding it into a complete flow. The qualitative 
behavior of solutions off the set of singularities is the same for both vector fields. The second ap- 
proach involves surgery and it is usually called regularization by surgery or block-regularization. 
Roughly, the idea is to excise a neighborhood of the singularity from the manifold on which the 
vector field is defined and then to identify appropriate points on the boundary of the region. 

Both approaches have been applied to the Kepler problem in Euclidean space. The second 
approach was first used for the Kepler problem by Easton ||, while the first approach has several 
variants. We mention only the most relevant for our work. As far as we know the regularization 



CN ' of the planar Kepler problem using the first approach was first discussed by Levi-Civita |T3|| 



Another beautiful incarnation of the second approach, due to Moser |jl7| , consist in showing 
that the flow of the n-dimensional Kepler problem on surface of constant negative energy is 
conjugate to the geodesic flow on the unit tangent bundle of S n . The main disadvantage of 
Levi-Civita and Moser's regularization methods is that they handle separately each energy level. 
This disadvantage is partially removed by a regularization procedure due to Ligon and Schaaf 



14|| . The Ligon-Schaaf regularization procedure allows to handle together all negative (resp., 
all positive) energy levels. However, negative, positive and zero energy levels still cannot be 
handled together with this procedure. The treatment of the regularization in the original article 
by Ligon and Schaaf requires laborious computations, and a somewhat simplified treatment of 
the Ligon-Schaaf regularization map is due to Cushman and Duisteermaat 0] and Cushman 
and Bates ||. Recently Marie |L5|, and Heckman and de Laat [Tl| posted on arXiv preprints 



that give another simplified treatment by showing that the Ligon-Schaaf map can be understood 
as an adaptation of the Moser regularization map. Regularization by surgery handles all the 
energy level at once, however this is a completely different approach that has other shortcomings. 
In fact, on the one hand it is difficult to find the attaching map and on the other hand this 
kind of regularization is not too helpful in understanding the global flow of the system and the 
near-collision orbits. 
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The Kepler problem on the three-sphere S 3 is the main topic of this paper. The Kepler 
problem and n-body problem on spaces of constant curvature are over a century old problems 
and recently have generated a good deal of scholarly interest. See |6|, [7|, [|, 19, 2C] for some recent 
results and some history. See [|J for a more exhaustive survey of recent results. 

The flow of the Kepler problem on the three-sphere S 3 is incomplete and several procedures 
have been used to regularize its vector field. The surgery approach to regularization was applied 



in |20 to the Kepler problem on a class of surfaces of revolution that include the two-sphere. 
The Levi-Civita and Moser's approaches were used in 



In the present article, we start by deriving the equations of motion using the method of 
Dirac brackets (section 2). In section 3, we describe the integrals of motion of the system and 
derive the equations of the orbit. In section 4, we describe the Moser regularization applied to 
the Kepler problem on S 3 . The details are given for the negative energy level sets, where the 
orbits lie completely in the upper hemisphere. The same approach applies to the part of any 
energy level set that lies in the upper hemisphere. Furthermore, we carry out the Ligon-Schaaf 
regularization directly (section 5) and compare it with the case of Euclidean Kepler problem. 
It turns out that the Ligon-Schaaf regularization for the negative energy part the two systems 
are naturally related by the gnomonic map (section 6). The gnomonic map takes the upper 
hemisphere to the full Euclidean space, and the two Kepler systems are related by a rescaling 
of the induced map on the tangent bundle. More precisely, let <3> and $ c be the Ligon-Schaaf 
maps for the Kepler problem on S" 3 and R 3 , respectively, and let \l/ denote the map induced by 
the gnomonic projection, then $ = $ c o ^. The map \1/ is not symplectic, which implies the 
non-symplecticness of the Ligon-Shaaf map for the spherical case. In the same section, we also 
show that the gnomonic transformation relates the Moser regularization of the two systems. 

The relations uncovered in this paper are summarized in the diagram below. 



Kepler on S 3 , 



9, Section [fj] 



Kepler on 




geodesic flow on S 3 



2. Preliminaries 

Let (•, •) be the Euclidean inner product in M 4 , and || ■ || the Euclidean norm in M 4 , and let •, 
x, and | • |, be the usual dot product, vector product and Euclidean norm in M 3 , respectively . 
Let q = (qo, q) with q = (qx, g 2 , (fe) and v = (v , v) with v = (vx, v 2 , V3) be canonical coordinates 

in TM 4 , with the symplectic 2-form to = dq. L A dvi. 

i 

Consider the following Hamiltonian system (H, TM 4 , ui) with Hamiltonian 



(2.1) 



H{q,v) = -(v,v) + V(q) 
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Qo 



where the potential energy V(q) is 

n,) = -7- (1 _ ?o2)1/2 

This Hamiltonian describes a particle that moves in M 4 under the influence of potential V(q). 

When the particle is constrained to move on the unit 3-sphere S 3 C M 4 , the system is restricted 
to the tangent bundle of the 3-sphere: 

TS 3 = {(q, v) e TR 4 \(q, q) - 1 = and (q, v) = 0} 

The Hamiltonian vector field Xh of H on TM 4 does not restrict to the Hamiltonian vector field 
Xh\ ts z of the constrained system. There are two approaches in computing Xh\ t$s using either 
a modified Hamiltonian function H* (resulting in Xh*, see Cushman ||) or a modified Poisson 
bracket, the Dirac-Poisson bracket {-, •}* (resulting in X H ) on TM 4 , such that the restriction of 
the resulting vector fields coincide with Xh\ ts3 - The two approaches are related since 

{F,G}*| TS a = {F*,G*}| TS 3 

The Dirac-Poisson bracket can be explicitly written down in this case. 

Lemma 2.1. The Dirac-Poisson structure {-,-}*\ts 3 is 
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or 



|T5 3 



TS 3 



Q3V2 ~ Q2V3 


UaiC[p}*\TS3 = 



where a, (3 = 0, 1, 2, 3 and L a p = qpv a — q a vp- 

Proof. The phase space TS 3 is given as a subset of TM 4 by the constraints 

ci(g, v) = (q, q) - 1 = and c 2 (q, v) = (q, v) = 
The Dirac-Poisson brackets are given by the relation 

{F, G}* = {F, G} + J2 C*{F, Ci }{G, Cj } 



1 -j 



where Cij are the elements of the inverse of the matrix with entries {q, Cj}. In this case 



C 



2(9, ?> 



The result follows from a computation. 

Lemma 2.2. The constrained equations on TS 3 are 

q = v 

v = -W(g) - ((v,v) - 



□ 



(2.3) 



W(<z))) q 
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Proof. We first compute the equations of motion for the Hamiltonian function H with respect 
to the Dirac-Poisson bracket 

q = {q,H}* 
v = {v,H}* 

then we restrict the resulting vector field to TS 3 . 

□ 



Let {60,61,62,63,64} be the standard basis for M 4 . Then for our choice of the potential, the 
equations of motion take the form: 



(2.4) 



q = v 

ie q 
v = - O, v) + 7- ji^J q 



(1 - qlfH 



(1"% 2 ) 3/2 ' 



3. Conserved quantities and the equations of the trajectory 

It is easy to see that the Hamiltonian H and the angular momentum fi = q x v are integrals 
of motion of the vector field Xh*\ts 3 = -^h\ ts3 - There is another interesting integral 

A Q 

A = 7T X \l — 7- — - 

m 

where 7r = q v — v q. The vector A is a spherical generalization of the vector known as the 
Runge-Lenz vector in Euclidean space. Priority, for the Runge-Lenz vector in Euclidean space, is 
sometimes attributed to Laplace, but appears to be due to Jakob Hermann and Johann Bernoulli 



lOfl . The only difference in form between this integral and the one for the Kepler problem in 



Euclidean space is that in the latter problem the velocity v appears in A instead of n. Note 
that, as in the case of the Kepler problem in Euclidean space, there is no universally accepted 
definition of the Runge-Lenz vector. The most common definition is given above, while the 

common alternative e = — is also called the eccentricity vector. 

7 

Proposition 3.1. The Runge-Lenz vector is an integral of motion. 



Proof. 



dA d-K q (q ■ q)q 

IT x M " 7fr + 7 1 , 3 

at dt \q\ \q\ 6 

, ■ ■ \ q (q ■ q)<i 

={q v + q v - qv - qv ) X /i - 7— + 7 



(3.1) 



v Q v - (v,v)q + ~f- 



?o 2 



13/2 



v lq-vq 

7|— T + 7 1 io , using 

\q\ 



q-v v-q 



7 / v Qo 



g| 3 /2 



X fj, 



1 



W2 W x 



7 



{v(q ■ q) - q{q ■ v)) 



7 



13/2 



[qx n-qx V) = Q 



□ 
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Remark 3.2. In analogy with the Euclidean case, there is an another conserved vector, the 
binormal vector, 



B = ^- i rr .a (a* x q) 
mm 



It can be shown that A = B x [i, and B is "binormal" because it is normal to both A and \i. 

Regard q, v, fi, A, B as vectors in M 3 . Since /x = q x v is a constant of motion, it points in a 
fixed direction. Moreover it is orthogonal to q and v. Therefore, when moving along an orbit, q 
varies on a plane (in M 3 ) orthogonal to \x. Furthermore, it can be shown that (A x B) ■ q = 0, 
A ■ /x = 0, and B ■ fi = 0. Hence, if A and B are both non-zero they form a basis for the plane 
where q lies (however, they are zero on circular orbits). 

Hence, each orbit must lie on a three-dimensional subspace. Moreover, since (q, q) = 1, the 
orbit is a curve on the two-sphere S 2 . If is used to denote the angle between q and the fixed 
direction of A then 

(3.2) A ■ q = \A\\q\ cos0 = q • (tt x /x) — 7 

\Q\ 

By permutation of the terms in the scalar triple product and note that n = q$v — v^q, 

q ■ (tt x n) = n ■ (q x 7r) = g M • (q x «) = ^0 1 A* | 2 
Let |q| = r, and go — -2- Rearranging equation (p.2|) yields 

COS ' 



(3.3) r V 7 

r 2 + z 2 = 1 

where (r, 0, 2) are cylindrical coordinates on the space spanned by A = (0, A), B = (0, B) and 
eo (where {eo, ei, e2, 63} is the standard basis for IR 4 ). 

Write the equation ( |3.3|) in spherical coordinates (p, 0, 6) using the formulas r = psin 9, = 0, 
and # = pcos#, we obtain 

r = 1 — p7 1 H COS 1 

tanb 1 \ 7 

This recovers the known formula for the orbits of the Kepler problem on the sphere |3| (compare 

also with [[0|). The orbits are always closed. Each orbit is obtained by intersecting a conical 

surface with vertex at the origin and a sphere, which reduces to a circle when \A\ = 0. In fact, 

the first of equations ( |3.3|) is a conical surface (and, in particular, one nappe of a quadric conical 
1-4-1 

surface for < 1), since for = O it reduces to the equation of a line through the origin, and 

7 

|_4I 

for z = Zq it reduces to a conic section. The constant e = lei = is called the eccentricity of 

7 

the orbit, which explains why e is called eccentricity vector. When e = 0, the orbits are circles, 
when < e < 1 they are curves corresponding to ellipses on the plane, e = 1 corresponding 
to parabolas, and e > 1 corresponding to hyperbolas. Note that the orbits are confined to the 
upper hemisphere when < e < 1, while it is not the case when e ^ 1. We refer the reader to 
H for a more detailed discussion of the orbits. 
The Hamiltonian H can be rewritten as 

H= l -{\^ + \rf) + V{q) 



6 SHENGDA HU AND MANUELE SANTOPRETE 

and the length of the eccentricity e vector is 

M 2 = l + M^(2#-|/x| 2 ) 



r 



There is another constant of motion 



which satisfies the following equation: 

(3.4) |e| 2 = l + 2^\n\ 2 

T 

Define the modified eccentricity vector e = —ue, where 
then the following hold 

/x ■ e = 
||/xf + ||£|| 2 = z/ 2 > 

where the second equation follows from (WM. These relations are equivalent to 



(/x + e) • (/i + e) = (/x - e) • (/x - e) = v 2 > 
which defines a smooth 4-dimensional manifold diffeomorphic to S 2 x S 2 .. 

4. HODOGRAPH AND MOSER'S REGULARIZATION 

In this section, we describe the Moser's regularization for the Kepler problem on S 3 . The 

approach follows closely Milnor's work for the Euclidean case. Without loss of generality, 

suppose that n = (0,0, fi). Then /i = q = ((ft, 0) and v = (vi,V2,0). Furthermore, 

suppose that A = (\A\,0, 0) 7^ and let <fi be the angle formed by q with A (as in the 

1^4] 

equation of the trajectory). Since A — B x /Lt, it follows that |-B| = -j — r(0, 1,0). Then /x x q = 



\fj,\ \q\(— sin0, cos <f), 0), and from the definition of -B, it follows 



7z = B + -P-r(u x q) = —4(0, 1,0) + (- sin 0, cos 0,0) = r^rf- sin 0, J — - + cos 0,0) 

Ia*I ImI ImI w 7 

I A\ 7 

Thus, for each orbit, n moves along a circle centered at - — -, with radius - — , and lying in the 



plane through the origin that is orthogonal to \x. 
Conversely, starting with a circle with the equation 



if , \ A \ 

7r = - — r — sin 0, h cos 0, 

m V 7 



a direct computation gives 



7 / A 
\x cos 6 = q fi = q x 77 = - — - sin 6 0, 0, 1 H cos 1 

Ia 1 ! V 7 
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Hence, we recover the equation of the orbit in spherical coordinates 



1 



tan 9 



_7_ 
1/4 



1^1 

7 



COS ' 



It follows that the hodocycle completely determines the corresponding Kepler orbit. We can 
now prove the following theorem. 

Theorem 4.1. Fixing some constant value of E = H — \n\ 2 /2 < 0, consider the space of 
all the vectors ir such that tt ■ 7r > 2E, together with a single improper point tv = oo. Such 
space possesses one and only one Riemannian metric ds 2 so that the arc-length parameter J ds 

/it 
— : — -. This metric is smooth 
Qo\q\ 

and complete, with constant curvature —2E, and its geodesies are precisely the circles or lines 
t — > n(t) associated with Kepler orbits. 

Proof. From (|4.1| ) it is clear that if E < the orbits must have go > and thus are limited 

to the upper hemisphere. In this case e = < 1 and only the space M% is relevant. In fact 

7 

and thus, if E < 0, | -rr 2 1 > 2E, since in such case go > 0. However, if 



E 



7T 



7" 



9o 



a - <? 2 ) 1/2 

E > things get more complicated. 

Using the equation of motions it is easy to show that 



7T = g t> + q Q v - v q - v q 



1Q_ 

kl 3 



and thus \n\ 



7 7 
- — -. Dividing this equation by the definition of the time rescaling ds = — ; — : dt 



and using the fact that 
(4.1) 

we obtain 



9b 9 



79o 

\q\ 

diz 



\*\ 2 -\{2H-\rf) 



ds 



(4.2) 



ds' 



= 7j9o| . 

\q\ 

4d7r ■ diz 



7T 



2 -E 

4d7r ■ diz 



(|tt| 2 - (2H- \ii\ 2 )) 2 (\n\ 2 -2E) 2 



Thus, there is one and only one Riemannian metric on the spaces which satisfies our 
condition, and it is given by formula (|4.2| ). 

To describe what happens in a neighborhood of infinity, we work with the inverted velocity 
coordinate w = (see, for example, Milnor's paper |16| for a discussion of inversion). Since 



the differential of w is dw = ^— L — we have dw ■ dw = — -. Hence 



(4.3) 



ds 



7T ■ 7T ) 

2 



7T • 7T 



4dw ■ dw 



(1 - 2Ew ■ w) 2 

The metrics given in equation ( |4.2|) and ( |4.3| ) have constant curvature —2E. 



□ 
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Since E < 0, the corresponding metric space is a round 3-sphere and t — > 7r(t) is thus a 
geodesic on a round 3-sphere. 

5. LlGON-SCHAAF REGULARIZATION 

5.1. so(4) momentum map. The components of the angular momentum and of an opportunely 
rescaled eccentricity vector form a Lie algebra under Poisson bracket which is isomorphic to 
50 (4). This gives the momentum map of the Kepler problem on the sphere. 

Let g = so (4). For a suitably chosen basis {X\, X2, X3, Y±, Y2, Y3}, the Lie bracket is 

(5.1) [X t , Xj] = [Y h Yd = e ijk X k and [X h Yj] = e ijk Y k 

The identification g = so (3) © so (3) can be seen via substitution -^{Xi + Yj) and -(X, — Yj). 

The basis can be thought of as coordinate functions on g*, then the Lie bracket defines a Poisson 
bracket {, } B on C°°(g*), which defines the Lie- Poisson structure on g*. 
We first write down the brackets of the components of fi and A: 

Lemma 5.1. The brackets of the components of fi and A are: 

(5.2) {/ij, fij}*\ TS 3 = eij k fj, k , {/ij, Aj}*\ TS 3 = eij k A k , {Ai, Aj}*\ TS 3 = — 2(H — \n\ 2 )eij k fj, k 

In the Appendix, we give a sketch of the proof of the Lemma above and of the Proposition 
below. 

Proposition 5.2. Let g = so (4), C(H) = H + + H 2 and 







2 + 


C{H) 




\A 


2 



Then the following is the momentum map of the Kepler problem on the sphere 

p=( t i,r ] (\ f x\ 2 ,H)A):TS 3 ^Q* 
It is a Poisson map with respect to the bracket {,}* on TS 3 and the Lie-Poisson bracket on g*. 

5.2. Delaunay Vector field. Let (x,y) be the coordinates on TM 4 = R 4 © M 4 . Then TS 3 C 
TM 4 is given by (x,y) = 0. Let T + S 3 = {(x,y) G TS 3 \y 7^ 0} be the tangent bundle of S 3 less 
its zero section and Cj = ulx+s 2 the restriction to T + S 3 of the standard symplectic form 00 on 
TM. 4 . Consider the Delaunay Hamiltonian on T + S 3 : 

(5 - 3) H -\if) 

which resembles the Kepler Hamiltonian written in Delaunay coordinates (see for example [0). 
It is clear that % is invariant under the standard action of 50(4) on T + S 3 . 
The integral curves of the Delaunay vector field X-^ satisfy 

dx 7 2 



(5.4) 



dt (y,y) 2 
dy 7 



y 



-X 



dt {y,y}' 

It can be proved that the Delaunay vector field Xy_ is a time rescaling of the Hamiltonian vector 
field of the geodesic flow on the unit sphere. The space so (4)* can be naturally identified with 
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/\ 2 (]R 4 )*. Under this identification, the momentum mapping of the standard action of SO (A) on 
T + S 3 is: 

J(x,y) = x Ay e /\ 2 (R 4 Y = so(4) 
A more detailed description of these facts can be found in |J. 

5.3. Regularization. All the orbits of the Kepler problem on the sphere that have E < can 
be regularized at once, with a slightly modified Ligon-Schaaf map. Let 

J : TS 3 ^so(4) : (q,v) \-> Qu,e) 

We begin our search for a Ligon-Schaaf map by noting that the image of J is the same as the 
image of J (see || for a detailed study of the Delaunay vector field and its momentum map). 
This suggests that the two maps are somewhat related, even though J is not a momentum map. 
Note that this situation differs from the case of the Kepler problem in IR 3 , studied by Cushman 
|J, where the maps that are related are the momentum maps. 

Theorem 5.3. The smooth map $ : (q, v) i— > (x, y) intertwines the momentum map J and the 
map J, that is <&*J = J, if an only if 

(5.5) $ : (q, v) \-> (x, y) = (a sin <p + (3 cos <p, v[—a cos </> + j3 sin ip)) 
where 

(5.6) a = (a , a) = ( q ■ 7r, — 7r 

\q\ 

7T ■ 7T 

7 



(5.7) /3 = (/3 ,/3) = 

tp is an arbitrary smooth real function and v 



Q_ 


q ■ 7r 


\q\ 


79o 




—A 


-i, 









Proof. Suppose $ intertwines the momentum map J and the map J. Write 
l 3 xl = IR 4 , y = (y , y)etfxl = R 4 . Then $*J7 = J is equivalent to 

x x y = q x v 

x y - y x = -y== ^7 - tt x (q x v)^J = Mq + Nir 

where 

1 / 7 7T • 7r\ , 1 / 77 ■ q 

(5.9) M = - and N 



V^2E \ \q\ q J V^2E \ q 

Suppose q x v = — (q x ir) ^ 0, then it follows that x and y lie on the same plane of q and 7r. 
Qo 

Since q and 7r are linearly independent, we obtain 

x = aq + bir 

(5.10) , 

y = cq + air 
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Since q x v = — (qx 77) = x x y = (ad — bc)(q x 77) we find ad — be = — . Substituting (|5.10|) 
into ( |5.8| ) and using the linear independence of q and 77 gives a set of linear equations for Xq 
and y . Since ad — be = — these equations may be solved to give 

x = (aN - bM)q 
y = (cN - dM)q 

Since (2, y) G T + S 3 , 

(5.11) l = x.x + x 2 

= x-y + x y 

Substituting the expressions for x and y and the expressions for xq and yo above into ( |5.11| ) 
yields 

, 2/ \ , ( T ' 9 . 7<?o ,\ 2 

(5.12) 1 = a (q ■ q) + a H 

2 



(5 ,3) = ac (, • , + + H + fc) ( 7 ^|) + M 

where we have used the identities 

(5-14) (,.*)- MiV^ = 7^|| 

77 • 7T + M q = 7 



|q| 2 (-2£) 

which follow from the definition of M and TV and the identity 

77 ■ 77 — 7- — r = 2.fc + 7- — - 

\q\ \q\ 

Multiplying ( p,12j ) by c and ( |5.13| ) by —a and adding the resulting equations gives 

_ 7(tt ■ q) _ b _ l 2 qo 



\q\(-2E) \q\*(-2E) 
Similarly multiplying (|5.12| ) by d and ( |5.13| ) by —a yields 



a = a q ■ q + -, r h 0— 

y >/ g |g|(-2£) 



All solutions of (|5.12|) are parametrized by 

1 

a = : 

\q\ 



smip 

(5.15) 



, \T2E\q\ * q . 
= cos if sin (p 

7<?o go7 
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where 9 is an arbitrary function. Substituting the previous equations in the expressions for c 
and d gives 

7 



COS 09 

\q\yf=2E 

, kl . . (w-g) 

a = — sm (p H -== cos ip 

qo q V-2E 



(5.16) 



Conversely, a calculation shows that a map $ of the form (|5.5|) intertwines the momentum 
map J and the map J. □ 

Corollary 5.4. $ intertwines E and the Delaunay Hamiltonian, that, is, $*"H = E. 
Proof. 

2{y,y) 2 7^] — a cos a? + p sm (p\ z 
since (a, a) = 1, (p 1 , /3) = 1, (a, /3) = 0, as it will be shown below. 

/„ «\ i o ^-g) 2 1 (7r-q) 2 (7r -tt) (7r-g) 2 
(a, a) = 1 - 2-j-j 1 5-5 1 Yl~ 

= 1 + 2 {7T ' q)2 f — - ^ V (-^(^-g) 2 
7 2 <?o V 2 l«| / 7 2 9o 



1 + 2 
1 



7 9o 7% v'qo 

I <gr| 2 (tt • 7r) / 7r • 77 7g (~2E) 



7 2 9o V ^ 191 



(a,/3}=-^ + i^(7r.g)(7r-7r) 



^0 ^79o 
7r-q |q| 



V • q)(7T • 77 ) =0 

□ 



The next result is a computation we will use to understand the relationship between the 
Kepler vector field Xh and the Delaunay vector field X^. 

Lemma 5.5. The derivatives of a and (3 along the flow generated by the vector field Xh are: 



da V^2E n d^ V^2E 
dt qo\q\ dt Qo\q\ 

Proof. Let t — > (q(t),v(t)) be an integral curve of the Kepler vector field Xh- Let a = 
a(q(t),v(t)) and (3 = f3(q(t),v(t)), where a = (ao, a) and (3 = (po,/3) are given by ( |5.6| ) and 
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( |5.7|) , respectively. Recall that 7r = q v — v q so that v — — (ir + v q). Moreover a simple corn- 
go 

diz *yq 

putation, shows that — — = —-, — nr. With these equations in mind, together with the expressions 

at \q\ 6 

for Xh, we compute 

da 1 1 

-dT = Wo {7 " q + 7 " v) ~^' q)qo 

'7 7T • 7T \ 1 1 

- r- r H H al 71 " ' <?) w o 2V 71 " ' 9)^o 

^go v l<?l go / v% ^% 

1 /(tt-tt) A 1 



and 



f?o|g| V <?o7 / ^gol<?l 



da i> (q ■ v) (tt • q) (tt ■ q) (n ■ q) . (it ■ q)7r 

— = — r-prQ 7T 7T 7T 2 V 

dt \q\ \q\ 6 79o IQo IQo IQo 

2"7T / 7T ■ 7T 7g " 



7?o V ' l 191 



- 71 " = — 



?olal v^go / qo\q 

Similarly, to verify the expression for -p we compute 

at 

d^o (w-ff)kl (7r-7r)|g| (w • 7r) (q • v) 

l 2 V0 + 



dt 7go 7go go7 \q\ 

(7T • q) (7T • 7r) (7T • 7r) , (7T • 71") 



go|<?l 2 7go 7gokl go7l<?l 

^ 7r • q ( 7r ■ 7r 7g \ _ \/—2E f \/—2E(7r ■ q) 



V 2 |q| / ?o|?| V 7go 



Qo\q\ 



-ot Q 



and 



d(3 \q\ir \q\7v (q ■ v)iz 

— T v o + _ 



dt vq vql vqo\q\ 



I T v H 2T — H 21 — ; 

7go 7go iQml i%\q\ 

y^2E ( q (q-Tz)iv\ V^2E 

-ex. 



go|g| 7go / go|<? 



□ 



We recall some standard terminology that we will use to describe the relation between the 
Kepler vector field and the Delaunay vector field. 
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Definition 5.1. Let X and Y be vector fields on the manifolds M and N, and let gf : M — > M 

and gj : N — >■ N be the corresponding flows. We say that X and Y (and the corresponding 
flows) are C k -conjugate if there is a C k -diffeoraorphism $ : M — > N such that 

We say that X and Y ( and the corresponding flows ) are C k -equivalent if there is a C ' -diffeomorphism 
$ : M — > N that maps the orbits of g x onto the orbits of g Y and preserves the direction of time. 
That is, there is a family of monotone increasing diffeomorphisms r x : R — > R such that 

Hg?(x))=gl {t) (<S>(x)) 
Remark 5.6. Note that differentiating Q(gf(x)) = 9^f x t)(^( x )) with respect to t yields 

(5.17) ^X(x) = f^j Y{9(x)) 

Hence, if this equation is satisfied, under the conditions of the definition above, the two vector 
fields are C k equivalent. 

The next theorem finally describes the relationship between the Kepler and the Delaunay 
vector fields. 

Theorem 5.7. The Kepler vector field Xh and the Delaunay vector field X% are C°° -equivalent 
with maps 

v) = (a sin ip + j3 cos ip, u(—a cos ip + j3 sin </?)) 

and 

dt 



to 



where ip = q ■ 7r(= a ). 

Proof. Differentiating 

(x(t),y(t)) = $(q(t),v(t)) = (a sin tp + f3 cos ip, v(—a cos <p + j3 sin <p)) 
with respect to t along a trajectory (q(t),v(t)) of Xh and using Lemma ^5| yields 

dx da dip d(3 d6 

— = — sin ip + a cos ip— — I — — cos w-psm ip— 
dt dt dt dt dt 

V^2E dip\ fo . . y/^2E fV^2E dip\ 
— (p sin ip — a cos <p) = I ; — ; — y 



and 



q \q\ dt J 7 V <?o|<?| dt 

7 d<p\ 



v Vz/^o \q\ dt J 



dy ( da d9 d(3 d6 

— = v \ — — cos B + asin ip— H — —smip + p cos 0— 
dt \ dt dt dt dt 



^2E dip 

q Q \q\ dt 

7 dip 

vqa\q\ dt 



(a sin <p + (3 cos ip) 



x 
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Taking <z> = (q ■ tz) it implies that 



dip 7 1 / v ■ 7T q ■ 7T q-7r.\ 7 



-go 



dt vqo\q\ v \ q q % J ^?o|<?| 

1/7T-7T q-7T 7 qr • 7r\ 7 

= - — t~ + — 2 rr ~ t; o — — i — i 

„2 



1 / O 79o\ 7 

2 I 7T ■ 7T 2— = - 



and thus by Corollary |5.4| we obtain 



*.{x H ) = I ^f y = \x H 



Let 



to «8(*) 

then ( [5.171 ) holds and hence the two vector field are smoothly equivalent. □ 

Note that the vector field of the Kepler problem in IR 3 is C fc -conjugate to the one of the 
Delaunay vector fields (Pj, 0]). And the diffeomorphism $ c intertwining the two vector fields is a 
symplectomorphism. However, the vector field Xh of the Kepler problem in S 3 is C fc -equivalent 
but not C fc -conjugate to the Delaunay vector field X H . 

6. Gnomonic Transformation 

The computations in the previous section are analogous to that of 0, |5|. However, there is 
a more direct link between the Kepler problem on S 3 and in Euclidean space. We first recall 
a classical result due to Appell and Serret @, ffifl , relating the Kepler problem on the upper 
hemisphere S* 3 to that on IR 3 via the gnomonic transformation. 

On the phase space ToIR 3 = (IR 3 — {0}) x IR 3 with coordinates (Q, V) and symplectic form 

3 

U3 = dQi A dVi consider the Kepler Hamiltonian 

i=i 

where ■ is the Euclidean inner product on IR 3 and \Q\ is the length of the vector Q. The integral 
curves of the Hamiltonian vector field Xh k on T 1R 3 satisfy the equations 

Q = V 

(6-1) T> Q 

V = ~ 1 W 

which (for 7 > 0) describes the motion of a particle of mass 1 about the origin under the 
influence of the Newtonian gravity. The momentum map of this system is 

Jk=(qxV,-^=(-Vx(QxV)--§- 
V V-2H K V7 \Q\ 
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This is a classical system and a detailed analysis can be found in ||, or in any good mechanics 
book. 

As before, S 3 is the unit sphere in M 4 . We now consider the gnomonic projection, which is 
the projection onto the tangent plane at the north pole from the center of the sphere: 

1 = Oo, q)^Q = 7= , 9 = — 

vi-|g| q ° 

The induced map on the tangent space at q is given by 

. v (q- v)q 1 qv 1 v 1 
v = Oo, v) i ^ + 3- = — v + ~^-q = —v - —q = — tt 

- \q 2 fi - lol 2 )2 % % qo q 2 % 



Theorem 6.1. The Kepler vector field Xu on S 3 (restricted to the upper hemisphere) and the 
Kepler vector field Xh k on IR 3 are C°° -equivalent with maps 

V:(q,v)^(Q,V)= (±ir 



qo 



and 

dt 



to 



q 2 (t) 

The map ^ together with t t— ?• r q is called the gnomonic transformation. 

Proof. Differentiating (Q(t), V(t)) = ^(q(t), v(t)) with respect to t along a trajectory (q(t), v(t)) 
of Xh and using Lemma [5.5| yields 

dQ 1 qv 1 

dt q q£ 

and 



dV _ dv \v\* (q-%) (q-vf 
-TT — qo-jT ~i 9 H 9 H 3 ( 

dt dt q q 



V l + 2 + 77; 2W2 + 9)) l"l 

<?o (i - %) ' J qo V 5o / 9o 

\q\ \Q\ qo\Qr 

where we used the fact that q^v 2 = (q ■ v) 2 . Consequently we obtain 



Let 




q 2 (t) 

then ( |5.17| ) holds and hence the two vector field are smoothly equivalent. □ 

The map \I/ of the theorem above has properties that are analogous to the map <£>. In fact 
it intertwines E and the Kepler Hamiltonian in Euclidean space, and it also intertwines the 
momentum map Jk of the Kepler problem in Euclidean space with the map J. We prove these 
properties below. 
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Proposition 6.2. Consider the smooth map \& : (q, v) — > (Q, V) 

(1) intertwines E and the Kepler Hamiltonian in Euclidean space, that, is, ^*Hk = E. 

(2) ^ intertwines the momentum map J% and the map J, that is ^*Jr — J ■ 

Proof. To prove the first part of the proposition we compute ^>*Hk- A simple computation 
yields 

(V*H K )(q,v) = \v-V-^-= l -{ql\v\* + 2 Q y + |<z|\ 2 ) 



since 



and 



2' |Q| 2— ' — u ' ™ V /T3^ 

o 77 ' 77 n 2 = 

2 V 1 - ?o 

v = g f H q = %v q = tt 

qo qo 

7 qo qo 

7^ = T 



\Q\ \q\ vi -qo 



fo 

We prove the second part by computing ^*Jk- A computation gives 

(V*J K )(q,v) = (q x V ,--^= (iv x(QxV)-j|JJ= (/*,e) = J(q,v) 

since 

Q X V = qr X + ^Tj^qJ = qxv = n 

IQI 7 V 1 - 9o T 

and ty*Hx = E by part (1) of the proposition. □ 

Recall that if (M, cum) and (N,un) are symplectic manifolds and / : M — > A is a diffeomor- 
phism then / is symplectic if and only if for all h, 

(6.2) f*X h = X hof 

see |]] for a proof. With this in mind we show that \& is not symplectic. 
Proposition 6.3. The map \1/ is not symplectic. 

Proof. To show that the map is not symplectic it is enough to find an Hamiltonian for which 
( |6.2|) is not satisfied. Let h = Hk, and let / = \&, then it can be shown that ( |6.2|) is equivalent 
to writing 

X Hlc {V(q,v)) = ^,X Hl<0 ^{q,v) = q*X E (q,v) 
The left-hand side of the equation is 

/ V \ / qov 

(6.3) fe(%«))= _o 

V 7 IQI a 

To compute the right hand side note that {fa, qj }*\ts 3 — e ijkqk and {fa, Vj}*\ TS 3 = e^Vk 
(see the Appendix), and thus {q, \h\ 2 }*\ts 3 — A* x q and {v, \n\ 2 }*\ TS s = /j, x v. Consequently 

Y ( q ,v) _ ( {q,H - |ji| 2 /2}*| TS s \ _ y(q,v) 1 ( MX q 




E " I KiJ-|/x| 2 /2}*| T5 3 ; "2 Ux« 
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where X^ ,v ^ denotes the part of the vector field corresponding to < q, H and 

. Differentiating Q(t) with respect to t along a trajectory of X E yields 



v,H ltd 



TS 3 



dQ v - \tx X q q[q ■ (v - \» X q)\ 1 _ x 1 , 

^ = Vl-kl 2 (X - |^P)3/ 2 =g(^ + (g-«)go 9) + ^ (-1*1 « + («■«)«) 

This gives the first three components of $,1^. Comparing with the first three components of 
X HK (V(q,v)) yields 

|q|2 _ i 

— a- (?ou + (<? ■ = Tj— (-|<7l 2,y + (9 • 

If i; is non-zero one can compare the coefficients of v. This comparison yields the impossibility 
1 = -, and hence the identity ( |6.2|) is not satisfied for h = Hk- □ 

Remark 6.4. The non-symplecticness of the map \1/ can be seen alternatively as follows. It can 
be shown that \1/ is the composition of the tangent lift with a scaling k : TS+ — > TS+ of the fiber 
direction by the factor of = 1 — |qr| 2 . As shown below, k is not symplectic with respect to the 
standard symplectic structure, which implies that the map *f> is not symplectic either: 

3 

K*dq A dv = dq A d((l - \q\ 2 )v) = (1 - \q\ 2 )dq A dv - 2 ^ v^dqi A c% ^ dq A dv 

Here, q is used as coordinates on 5 3 ; via the projection along the q direction. 

6.1. Relation to the Ligon-Schaaf regularization for the Kepler problem in M 3 . Recall 
the description of the Ligon-Schaaf regularization for the Kepler problem on IR 3 as given in f|. |5|] . 
The symplectomorphism $ c intertwining the vector fields of the Kepler problem and one of the 
Delaunay vector fields is given by 

$ c : (Q, V) i y (x c , y c ) = (a c sin y? c + C cos ip c , v c {-a c cos (f c + (5 C sin ip c )) 
7 

where v c = , y? c = a Cj0 and 

sJ—lti K 

(*c = ("c,0, ol c ) = —V Q, — V 

\Vc \Q\ 7 

ft = (/3 Ci0 ,/3 c ) = (^(V ■ V) - 1, i^iv 
V 7 

This calculation proves the following 
Proposition 6.5. $ c o = $ 

A straightforward consequence of the proposition above is the following 
Corollary 6.6. $ is not symplectic. 

Proof. Suppose $ is symplectic. Then we can rewrite the proposition above as \& = o $. 
Since and $ are symplectic it follows that \l/ is symplectic. This is in contradiction with 
Proposition □ 
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6.2. Relation to Moser's regularization for the Kepler problem in M 3 . First we recall 
Theorem 2 in jTj| for the Kepler problem in IR 3 . For a given constant value of h = Hr, define 

the space M h = {V\V -V > 2/i}U{oo}. Define also a Riemannian metric ds 2 = -j— — — — 
on M h . The arc-length parameter j ds of this metrie along any veloeity circle r -» V(r) is 
equal to the parameter J — -, where r denotes the time. 

Proposition 6.7. The Moser's regularization for the Kepler problems in M. 3 and the upper 
hemisphere are related by the gnomonic transformation. 

Proof. The gnomonic transformation is given by (Q, V) = ^(q, v) = ( — , n J and dr = — along 

the trajectory It's clear that \I/ maps to and the metrics correspond. The arc- 

f dr f § f dt 
length parameter becomes / -—=-7 = / — — = / — j — r. □ 

J \Q\ J SL J qo\q\ 

go 



Appendix 

Proof of Lemma \5.J\ . We verify only the third equation in ( p.2| ). Using Lemma |2.1| we obtain 

{A*i) QjYWs 3 — ^ijkQk, {HiiVjYWs 3 — CijkVk, {^i , }* \ TS 3 = e ijk^k 

and 



{7Ti,7rj}*| T s3 = qtVj - qjVi, \q\ Y\ts 3 = 



1 

\q\ 



TS 3 



Using the bilinearity property of Poisson brackets, expand 



to obtain 



It helps to recall the identity 



{Ai,Aj}*\ TS 3 = -2(H - \n\ 2 )e ijk fi k 



^■ijk^ilm OkmOjl ^jm^kl 



\TS 3 



□ 



Proof of Proposition |<5.ij . Let E = r](\fi\ 2 , H)A, and let E = (Ei, E 2 , E 3 ), then the proposition 



amounts to showing that {fix, /i 2 , A*3, E\, E 2 , E 3 } satisfies (|5.1|) with the Poisson bracket {, }*. 
We start with ( |5.2|) and compute the brackets of the components of E. For r] = r/(|/x| 2 , H), 

{iqA u iqA ] Y\ TS z=r l 2 {A i ,A ] Y\ TS z+r 1 {A i ,r 1 Y\ TS zA ] +T 1 {r 1 ,A j Y\ Ts3 A l 
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r a i*i . 9tj drj d\n\ 2 drj 

{Vi A/ \ts 3 — £ijkA k — — Oji — eu k A k —, — rr— — — —lei k i[iiA k - 



dfii d\n\ 2 d\fji\ 



Consequently 



on 

{rjAi^AjYlrss = ^{A,, Aj}*\ TS a + 2r]—^ ((// x A) t Aj - (/x x A)jAi) 

o\ii\ 



Simple computations show that 

n 2 (H - I if.l 2 W- ,.//.,_ - 



{r}Ai,r}Aj}*\TS3 = -2rj\H - \n\ 2 )e ijk n k - -^—^\ A\ 2 e ijk ^ k 



-2r ] 2 (H-\„\ 2 )--^ 2 ( 1 2 + 2H\„\ 2 -(\„\ 2 ) 2 ) 



Since \ A\ 2 = 7 2 + \n\ 2 (2H - |/z| 2 ), we have that |gL = 2{H - |/i| 2 ). Therefore 

~ ~ d\tiA\ 2 d 

(6.4) {-^i, £y}*|r53 = {77A, 774j}*| t <j3 = — e ijfc ^ fc = -(e ijfc /i fc )-— -|e| 2 

C|/i.| c I A* I 

Lastly, we determine the function 77. The equation (|6.4|) implies that 77 has to satisfy 

d\yA\ 2 

d\v\ 2 

Integrating we obtain 7/ 2 |A| 2 = — \fi\ 2 + C(H), and if \A\ 2 7^ we have 

v 2 

where C(H) is an arbitrary function of H. For \A\ 2 = we have \fi\ 2 = C(H). Substituting in 
the expression \A\ 2 = r y 2 + \^i\ 2 (2H — \n\ 2 ) yields 

7 2 + C(H)(2H -C{H)) = 







2 + C(if) 






2 



then the positive solution gives C(H) = H + y 7 2 + -f^ 2 - □ 
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